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We present a new contribution of the R-parity violating {^p) supersymmetry (SUSY) to neutri- 
noless double beta decay (0!//3/3) via the pion exchange between decaying neutrons. 

The pion coupling to the final state electrons is induced by the I^p SUSY interactions. We have 
found this pion-exchange mechanism to dominate over the conventional two-nucleon one. The latter 
' corresponds to direct interaction between quarks from two decaying neutrons without any light 

hadronic mediator like vr-meson. 

The constraints on the certain SUSY parameters are extracted from the current experimental 
^ I Oi^/3/3-decay half-life limit. These constraints are significantly stronger than those previously known 

or expected from the ongoing accelerator experiments. 
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Neutrinoless double beta decay (Oi'/?/?) has long been recognized as a sensitive probe of the new physics beyond 
the standard model (SM) (see [Q- j|]). Various mechanisms of OvPf] decay were proposed and studied in the last 
^ two decades. The conventional mechanism is based on the exchange of a massive Majorana neutrino between the two 
decaying neutrons. A new mechanism was found within supersymmetric (SUSY) models with R-parity violation (-^p) 
, in 0. {Rp = (—1)3^+^+25' -ftrhere S, B and L are spin, baryon and lepton numbers.) It was later studied in more 

■ details in Q. A complete analysis of this mechanism within the minimal supersymmetric standard model (MSSM) 
^—1 ' was carried out in jsj. 

\^ The nuclear Oi'/J/^-decay is triggered by the Qvj3[3 quark transition d + d u + u + 2e^ which is induced by certain 

0^ fundamental interactions. It was a common practice to put the initial d-quarks separately inside the two initial 
(—! • neutrons of a Oz^/3/3-decaying nucleus. This is the so called two-nucleon mode of the Oj^/3/3-decay [see Fig. 1(a)]. If 
O 1' the above OvPP quark transition proceeds at short distances, as in the case of J^p SUSY interactions, then the basic 
' nucleon transition amplitude n + n—fp + p + 2e~ is strongly suppressed for relative distances larger than the mean 
nucleon radius. 

\_ In this letter we propose a new pion-exchange SUSY mechanism which is based on the double-pion exchange 

■ between the decaying neutrons [Fig. |l|(b)]. At the quark level this mechanism implies the same short-distance I^p 

', MSSM interactions as in [|j. However, it essentially differs from the previous consideration of the SUSY contribution 
to the Oz^/3/3-decay at the stage of the hadronization. We assume that the J^p MSSM quark interactions induce 
$-H ' TTTT — > 2e transition at the middle point of the diagram in Fig. |l|(b). The importance of the pion-exchange currents in 
5^ 1 Oi//3/3-decay was first been pointed out by B. Pontecorvo Latter, this idea was quantitatively realized in 0- |^ 
for the case of the heavy Majorana neutrino exchange. It was shown that the pion-exchange contribution can not be 
neglected in this case. We will show that in the case of the I^p MSSM induced quark transition the pion-exchange 
contribution absolutely dominates over the conventional two-nucleon mode. 

The i?p-violating part of the superpotential breaking lepton number conservation is 

W^^p = KjkLiLjEk + X'ijf.LiQjDk. (1) 

Here L, Q are lepton and quark doublets while E, D are lepton and down quark singlet superfields. Indices i,j,k 
denote generations and Xijk = —Xjik- In what follows we concentrate on the so called "direct" SUSY contribution 
to the OvPP 1^- [|| depending only on the A' term of the superpotential in Eq. (^. The combination of both A' and 
A terms may lead to the "indirect" SUSY contribution accompanied by the neutrino exchange H], which wc do not 
consider in the present letter. 

Starting from the A' term in Eq. (^ the following effective quark-electron vertex has been derived [^ : 

=S^e-(l+,.)e^ 



(2) 



1 



These interactions violate the electron number ALe = 2. They are induced by heavy SUSY particles in a virtual 
intermediate state. An example of the Feynman diagram contributing to Cqe is given in Fig. ^ The color-singlet 
hadronic currents in Eq. (^) are Jp — u'^j^da, Js — u°'da, Jj^ — u°'a'^'^{l + j5)da, where a is a color index. 
The lepton number violating parameters rj in Eq. (|^) can be written in the following form 

\ rrig A J 
rn„ 3 / rriq 



^ AM2a,^ + -a2^^ C . (4) 

J \ Tin- / m . \ m- I I 



Here A = {V2TT/3)X'mGp^m- ^ and C = &{mq/meYe\^ - tm^Le - ^Lu^Rd{me/mqY - ^Lu^Lc- ^2 = gl/ (47r) and = 
f/|/(47r) are SU{2)l and SU{3)c gauge coupling constants; and are masses of the gluino g and of the lightest 
neutralino x- The latter is a linear combination of the gaugino and higgsino fields x = ct^B + (i^W^ + S^H^ + 7j(.-^2 ■ 
Here and B are neutral SU{2)l and U gauginos while ff^ ? -^^i are higgsinos which are the superpartners of the 
two neutral Higgs boson fields and H2 with a weak hypercharge Y = —1, +1, respectively. The mixing coefficients 
Uj^T I3^,^xt^x '^^'^ be obtained from diagonalization of the 4x4 neutralino mass matrix. Neutralino couplings are 
defined as tLtp — — '?3('(/')/3x + tan^vF (Tsiip) — Q{ip)) a^, tR^f, — Q(V') t&nOwctx [@- Here Q and T3 are the electric 
charge and weak isospin of the fields ip — u,d,e. In Eqs. we used the universal squark mass rriq ansatz at 

the weak scale ruu ~ ?7ij « niq. This approximation is justified by the constraints from the flavor changing neutral 
currents and is sufficient for our analysis. 

Now we have to reformulate the quark-lepton interactions in Eq. (|^) in terms of the effective hadron-lepton interac- 
tions, which is necessary for the further nuclear structure calculations. The effective Lagrangian taking into account 
both the nucleon (p, n) and 7r-meson degrees of freedoms in a nucleus can be written as follows: 

£- pF^n pr«n e(l + 75)e'= 



2 m 



G 



F ™4 



m^a^ (tt Y e(l -1-75)6'^ 



2™ 

+9s P T^'^ ■ (5) 

Here C2Ni t^-Ke and C^^f^ describe the conventional two- nucleon mode, pion-exchange mode and pion- nucleon interac- 
tions, respectively. They correspond to the first, the second and third terms of the second part of the equation (||). 

= 13.4 ± 1 is known from experiment. The two-nuclcon mode contributions C2N to the 0^/3/3-decay with different 
operator structures F^*) were derived and studied in [Q- [|| within the I^p MSSM. 

In this note we concentrate on the effect of the pion-exchange term C-^e- The basic parameter a,r of the Lagrangian 
£,re can be approximately related to the parameters of the fundamental Lagrangian Cqe using the on-mass-shell 
"matching condition" < 7r+|£ge|7r^ >=< Tr+l^jreK" >• The solution of this equation is a^r = ji'Uq + Vf){C'p + ^s) — 
■^TjqC^, where < Tr~^\JiJi\TT~ >= — m^Ci with i = P,S,T. Thus, we obtain the approximate hadronic "image" Cyre of 
the fundamental quark-lepton Lagrangian Cqe given in Eq. (||) . 

The contribution of the Jp,s,T currents to a,r can be estimated within the vacuum insertion approximation (VIA). 
Applying Partial Conservation of Axial Current (PCAC) we obtain 

g 

< TT+I J_pJ_p|7r" > =: - < TT+I Jp|0 >< 0| Jp|7r" > 
o 

= -^L-f xI = -™,rCp- 6 

3 [mu + mdy 

where 8/3 is a combinatorial color factor and /^r — 0.668 m^. Taking the conventional values of the current quark 
masses niu — 4.2 MeV, nid — 7.4 MeV one gets Cp « 342. Within the VIA we have c,, = = since < 
0| J5|7r(p^) >=< 0| J^''|7r(p7r) 0. The scalar current matrix element vanishes due to the parity arguments, the 
tensor one vanishes due to J^"" — —Jt^ and the impossibility of constructing antisymmetric object having only one 
4- vector p^. Thus, we expect the Jp contribution to be dominant. 

The Jp dominance also follows from the non-rclativistic quark model (QM) Within this model one can 

calculate < 7r+|JpJp|7r^ > using the closure approximation for the intermediate meson states W. After quite 
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tedious calculations we end up again with Cp >> c^^. In this case the numerical value Cp « 1100 is larger than 
that in Eq. (^) since in addition to the vacuum state there are other intermediate states taken into account. In what 
follows we use both the VIA and the QM values of Cp . 

The large coefficient Cp enhances the pion-exchange contribution to the Oi^/?/?-decay. This enhancement factor is 
a generic property of the I^p SUSY- models generating at low energies the JpJp interactions [see Eq. (||)]. There is 
another factor enhancing the pion-exchange contribution compared to the two-nucleon mode. As explained latter on, 
it stems from the fact that the pion-exchange is longer ranged and thus covers a larger interval of the intcr-nuclcon 
distances enhancing the nuclear matrix elements over the two-nucleon mode. 

Starting from the Lagrangian Che in Eq. (|^) it is straightforward to calculate the contribution to the Qvl3(3- 
matrix element TZovpfi which corresponds to the fundamental vertex Cq^ in Eq. (H). It consists of the two terms 
T^Oupp — ^0^3/3 + '^Zvpp describing the conventional two-nucleon mode Ti-^^pp and the pion-exchange contribution 
'^OvBB- The relevant Feynman diagrams are given in the Fig. |l|. The corresponding half- life formula reads 



Mf + 7?; MY -f (r^q + n^) M^^ ' . (7) 

Here Gqi is the standard phase space factor tabulated for various nuclei in 1^ and iriA = 850 MeV. The analytic 
form of the two-nucleon mode nuclear matrix elements given in Here we present the new pion-exchange 

nuclear matrix element defined as 

JTL 

M^"" = {Mgt,. + Mt,.) , (8) 

rrie 

The partial Gamow- Teller and tensor matrix elements are 

Mgt,. = < 5] r+T+a,, (—) F,{x^)\0+ >, (9) 



where 



Mt,. = < 0+1 ^r+r+ 5,, (—) F,{x.)\Ot >, 



(10) 



Sij — 3(7j • T-ij (7j • Tij (Ji • fJ J , (J{j — (7 i • <7j^ 

nj = {n - rj)/\ri - fjl, n.j = \ri - rj\ (11) 

and — m-j^rij. Here fi is the coordinate of the "i-th" nucleon. The pion-nucleon structure coefficient in Eq. (^) is 
given by 

1/ \2/ \4/ x2 



96 V "^p / V "^A / \JaJ 
where Ja — 1.261. The pion-exchange SUSY potentials are 

Fi(a;) = (x-2)e-^ F^ix) ^ {x + l)e-\ (13) 

The most stringent experimental lower limit on the Oi^/3/3-decay half- life has been obtained for ^^Ge |l3); what 
favors especially this nucleus for nuclear structure calculations. In this letter we pay our attention only to this 
isotope. We have employed the renormalized quasiparticle random phase approximation with proton-neutron pairing 



(fuU-RQRPA) 1 12 1 to calculate both the two-nucleon and pion-exchange nuclear matrix elements governing the I^p 
SUSY Oi^/3/3-decay of ^^Ge. The fuU-RQRPA includes the Pauli effect of fermion pairs and does not collapse for a 
physical value of the nuclear force strength. To include the Pauli principle more correctly we do not use the quasi- 
boson-approximation to derive the Quasiparticle Random Phase Approximation (QRPA). If one includes the exact 
Fermion commutation relations for nucleon pairs (two quasiparticles) as a QRPA expectation value, one obtains the 
renormalized QRPA (RQRPA), which is stable against the collapse of 2v[3l3 Gamow- Teller transition. Therefore the 
RQRPA offers a significantly more reliable treatment of the nuclear many-body problem for the description of the 
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Oi//3/3 decay. Thus it also allows one to establish more reliable constraints on the SUSY parameters from the best 
available experimental lower bound on the Oi^/3/3-decay half-life. We have found the following numerical values of the 
nuclear matrix elements for ™Ge: Mf^ = -61; My = 0.85; M'"^ = -1800(QM), -600(VIA). The pion-exchange 
matrix element is given for QM (non-relativistic Quark Model) and VIA (Vacuum Insertion Approximation) values 
of the coefficient Cp . It is apparent that in both QM and VIA cases the dominant contribution to Eq. comes from 
the pion-exchange mechanism corresponding to M.'^^ . The VIA value M.'^'^ = —600 we will use for conservative 
estimations. It is worthwhile to note that the above nuclear matrix elements are quite stable with respect to variation 
of the nuclear model parameters. The uncertainty of the calculated values of -^^^ a-nd M"^^ does not exceed 20%. 

Now we are ready to extract the constraints on the I^p MSSM parameters from the non-observation of Oz^/3/3-decay. 
The current experimental lower bound on the '^^Ge Oz//3/3-decay half-life is T^'j^'^"'''''^ {{)+ 0+) > 9.1 x lO^^ 

years 90% c.l. Combining this bound with Eq. (Q) and the above given numerical values of M'^^ we get a constraint 
on the sum of the effective MSSM parameters rjg + rjj < 2.1 x 10~^. If one does not include the pion-exchange 

contribution then one gets a constraint rjq + 0.014?;^ < 7.8 x 10^^ from the remaining 2N-mode. It is essentially less 
stringent than the above given 7rN-mode constraint by more than one order of magnitude for rjg and by three orders 
for 77y. 

The gluino and neutralino contributions to rji can not cancel each other within the present experimental limits on 
their masses and couplings [see j^]. Therefore, we can extract from the above limit on rji the constraints on these 
individual contributions. The gluino contribution constraint is 

1/2 



a;„ < 2^0(1,18) X lO-(^) (^) , (14) 



for the VIA(QM) value of the pion matrix element parameter Cp . The neutralino contribution constraint is more 
complex because it involves more parameters: neutralino mixing coefficients, selectron and squark masses. However, 
it can be cast into the form of Eq. (jlj) under the phenomenologically viable simplifying assumptions. Assume that 
the neutralino is B-ino dominant a^^ >> jS^^S^^j^ and that rriq > nie ■ Then we get 



a;„< 5.2(3.07) xl„-.(^) = (^)"=. (15) 



If all SUSY particle masses in Eqs. (14)-(|15[) were at their present experimental lower bounds ||1J| rriq > 90 GeV, 
me > 45 GeV, mg > 100 GeV, > 19 GeV, we could estimate the size of the l^p coupling constant X'm < 
4.6(2.7) X 10~^. A conservative bound can be obtained using the SUSY "naturalness" upper bound mg q_^ < ITeV. 
It gives Ajii < 6.3(3.7) X 10~2. This limit and those in Eqs. (|l|)-(|l|) are the best known limits on the l^p coupling 
X'lii [see y and references therein]. 

It is interesting to compare the 0:^/3/3-decay and accelerator experiments from the point of view of their sensitivity 
to the l^p SUSY signal. Previously [^ it was shown that even in the 2N-mode the constraints from the Oi^/3/3-decay 
exclude the domain of the ^p MSSM parameter space accessible for the ongoing experiments with the ZEUS detector 
at HERA. This conclusion touched upon the region of the so called I^p resonant single squark q production mechanism 
in deep inelastic ep-scattering [ p^ . Taking into account the 7rN-mode makes this conclusion much stronger. The 
reason is that the excluded region becomes significantly larger than in case of the 2N-mode alone. Now this region 
extends so far that it would include the HERA domain even if the experimental lower bound on T^J.^^ was by about 
(conservatively) five orders of magnitude less than the currently existing one. 

Nevertheless, there is still a window for the HERA experiments in the region corresponding to another mechanism 
assuming i?p-conserving e + q production and their subsequent I^p cascade decays. This region is + rriq < 205 
GeV and Ajjj > 10^^ [QQ- the present letter we would like to stress that the 0j//9/3-decay constraints given in 
Eqs. (0)-(|l5[) dramatically reduces the above quoted region. Only a narrow part of it corresponding to very low 
values of the l^p coupling constant 10^^ < A'^ < 0.76(2.4) • 10~'^( ^™Q"^y )^ and mg -f- rriq < 205 GeV remains not 
excluded by the Oz^/?/3-decay constraints. Here we used the conservative VIA value of the pion matrix element cp and 
put m^=mg = 100 GeV (1 TeV). 

Summarizing, we point out that the SUSY contribution to the Oz//3/3-decay comes dominantly via the pion-exchange 
mechanism considered in the present letter. The conventional two-nucleon mechanism corresponding to 

nn — > ppee transition without light particles (pion or neutrino) in the intermediate state, brings only a subdominant 



^Note that unlike Oi^/S/S-decay searches HERA can probe within this mechanism I^p coupling constants other than A'm. 
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SUSY contribution. On practice the pion-exchange mechanism considerably enhances the sensitivity of the Qi'f3(3- 
decay to the supersymmetry. This aUowed us to obtain presently the most stringent limitations of the certain first 
generation I^p MSSM parameters. 
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